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Abstract 

We  study  the  limiting  behaviour  of  large  systems  of  two  types  of  Brownian  particles 
undergoing  bisexual  branching.  Particles  of  each  type  generate  individuals  of  both  types, 
and  the  respective  branching  law  is  asymptotically  critical  for  the  two-dimensional  system, 
while  being  subcritical  for  each  individual  population. 

The  main  result  of  the  paper  is  that  the  limiting  behaviour  of  suitably  scaled  sums  and 
differences  of  the  two  populations  is  given  by  a  pair  of  measure  and  distribution  valued 
processes  which,  together,  determine  the  limit  behaviours  of  the  individual  populations. 

Our  proofs  are  based  on  the  martingale  problem  approach  to  general  state  space  processes. 
The  fact  that  our  limit  involves  both  measure  and  distribution  valued  processes  requires  the 
development  of  some  new  methodologies  of  independent  interest. 
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Introduction 


In  this  paper  we  shall  study  the  limiting  behaviour  of  large  systems  of  two  populations  of 
Brownian  particles  undergoing  bisexual  branching.  As  an  example,  consider  the  following 
system: 

Assume  that,  at  time  t  =  0,  n  >  1  particles  of  two  types  (type  1  and  type  2)  are  located 
in  Rd  and  begin  moving  as  independent  Brownian  motions.  Each  particle  of  type  i,  i  =  1,2, 
dies,  independently  of  the  others,  after  an  exponential  time  with  mean  1/n.  Assume  that, 
at  the  time  of  death,  each  particle  is  replaced,  on  average,  by  one  particle  of  each  type, 
with  overall  probability  and,  with  probability  |,  by  nothing.  The  replacement  particles,  if 
there  are  any,  perform  independent  Brownian  motions  starting  at  the  point  of  death  of  their 
parents,  and  this  story  of  diffusion  and  (critical)  branching  continues  until,  with  probability 
one,  there  are  no  particles  left. 

Let  z\n\t ,  •)  be  the  measure  valued  processes  describing  the  positions  of  the  particles  at 
time  t.  That  is,  let 

Z-n\l,A)  =  Number  of  particles  of  type  i  in  A  at  time  t. 

Our  interest  lies  in  establishing  the  joint  behaviour  of  the  Z\n\t )  in  the  infinite  density, 
n  — >  oo,  limit. 

This  problem,  in  the  absence  of  particle  motion,  has  already  been  studied  by  Kurtz  [5]. 
The  addition  of  the  particle  motion,  however,  makes  the  problem  considerably  more  complex, 
and  rather  interesting. 

In  order  to  describe  a  typical  result,  define  the  following  two,  rescaled,  processes. 

*(.)W  =  +  yW(()  _  zj-'to  - 

n  n 

We  shall  study  the  limiting  behaviour  of  and  as  n  — >  oo.  This,  clearly,  will  tell 
us  about  the  limiting  behaviour  of  the  Z-n\  What  we  shall  find  is  that  while  A'*"*  has  an 
interesting  nontrivial  limit,  converges  to  zero.  That  is,  in  the  infinite  density  limit  the 
proportions  of  particles  of  each  type  are  identical. 

This,  naturally,  leads  one  to  try  to  re- rescale  the  difference  Y^n\  so  as  to  obtain  a  fluc¬ 
tuation  result  describing  the  rate  at  which  balance  between  the  two  populations  is  achieved. 
It  turns  out  that  this  is  best  done  via  a  third  process,  defined  by 

W{n)(t)  =  y(">(t)  +  t'nYln){s)ds. 

Jo 

We  shall  show  that  has  a  nice  limit,  as  n  — ►  oo,  as  a  (Schwartz)  distribution  valued 
process,  and  that  the  convergence  is  joint  with  that  of  A'*nL 

The  main  result  of  the  paper,  Theorem  3.3,  gives  the  details  of  this  convergence,  under 
a  more  general  setup  than  that  just  described. 

Since  our  study  involves  both  measure  and  distribution- valued  processes  we  shall,  un¬ 
fortunately,  need  to  start  with  some  technical  results  about  weak  convergence  for  measure 
cross  distribution  valued  processes.  Some  of  these  results  should  be  of  independent  interest. 
This  is  Section  1  of  the  paper. 
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In  Section  2  we  begin  describing  our  system  in  detail,  and  describe  also  the  previous 
result  of  Kurtz  [5]  noted  above.  Section  3  contains  the  main  result  of  the  paper,  and  Section 
4  is  devoted  to  proofs.  There  is  an  appendix  containing  the  proof  of  the  existence  and 
uniqueness  of  the  solution  a  particular  non-linear  evolution  equation  which  is  used  in  the 
proofs  of  Section  4. 

1  Preliminaries,  weak  convergence 

In  this  section  we  give  a  brief  introduction  to  the  martingale  problem  approach  to  general 
state  space  process.  We  start  with  some  notation  and  definitions.  Let  E  be  topological  space, 
B(E)  (C  ( E ))  be  the  set  of  bounded  (bounded  continuous)  Borel  measurable  functions  on  E. 
Denote  by  B(E )  the  <7-algebra  of  Borel  subsets  of  E  and  by  V  (E)  the  set  of  Borel  probability 
measures  on  E.  Let  A  be  subset  (not  necessary  linear)  of  B(E)  x  B(E). 

Definition  1.1  By  a  solution  of  the  martingale  problem  for  A  we  mean  a  measurable  stochas¬ 
tic  process  X  with  values  in  E  defined  on  some  probability  space  (Q,IF,P),  such  that  for  each 
(fid)  €  A  the  process 

(1.1)  /(*(*))-  fg{X{s))ds 

Jo 

is  a  martingale  with  respect  to  the  filtration 

{X(u))  du  :  s<t,  he  £(£)) 

where  Tfi  =  d(X(s)  :  s  <  t). 

Note  that  if  X  is  a  right  continuous  process  then  Tfi  —  Tfi . 

Definition  1.2  When  an  initial  distribution  p  e  V  {E)  is  specified,  we  say  that  a  solu¬ 
tion  of  the  martingale  problem  for  A  is  a  solution  of  the  martingale  problem  for  {A,p)  if 
PX{  0)-1  =  p. 

Let  be  a  complete  filtration.  Let  C  be  the  space  of  progressive  (i.e.  (.^-progressive) 
processes  Y  such  that  supt>0  £[|K(f)|]  <  oo.  Set 

(1.2)  IM  =sup£[|y(()|]. 

(>0 

We  do  not  distinguish  between  C  and  the  Banach  space  of  equivalence  classes  in  £,  deter¬ 
mined  by  the  norm  (1.2).  (  X  ~  Y  if  ||^  —  K||  =  0  ).  We  define  the  semigroup  of  operators 
J{s)  on  £  by 

(1.3)  (jwr)«  =  £[r<‘ +  .)!*]• 

Define 

A  =  l(Y,  Z)  e  £  x  £  :  K(t)  —  /  Z(s)ds  is  .Fj-martingale  j  . 


(1.4) 


A  is  called  the  full  generator  of  J{s). 

Let  G(,De[0,  oo))  be  the  cr-algebra  generated  by  the  simple  cylindrical  subsets  of  Z)£[0,  oo); 
i.e. 

G(De[  0,  oo))  =  a  (K itmy'  ( B(E)fm €  [0,oo),  metf), 

where  for  points  ti,t2,.  .  .,fm  in  R+,  the  projection  ntl<t2 . tm:  Z)£[0,  oo)  — ►  Em  is  defined 

by 

. tm(X )  =  (l(<l),*(<2),...,x(«m)),  V*  G  ££[0,Oo). 

When 

(1-5)  S(D£[0,oo))  =  S(££[0,oo)) 

we  shall  say  that  E  has  equivalent  Borel  and  cylindrical  cr-algebras. 

We  now  give  an  analogue  of  the  Theorem  4.8.10  [5]  without  the  assumption  that  the 
space  E  is  metric.  The  result  will  be  crucial  for  our  later  needs,  as  our  main  result  will, 
ultimately,  follow  directly  by  checking  that  the  conditions  of  this  result  are  satisfied. 

Theorem  1.3  Let  ( E,r )  be  a  completely  regular  topological  space  with  equivalent  Borel  and 
cylinder  o-algebras  on  Z)£[0,oo).  For  each  n  >  1  let  {Ef}  be  a  complete  filtration,  and 
let  A  C  C  ( E )  x  C  (E)  and  v  G  V  ( E ).  Suppose  X„,  n  =  1,2, . . is  a  Ef -adapted  process 
with  sample  paths  in  Z)£[0,  oo),  {A„}  is  relatively  compact,  and  PXn( 0)_1  =S>-  i/,  as  n  — ►  oo. 
Suppose  furthermore  that  for  each  ( f,h )  G  A  and  T  >  0,  there  exist  (£n,  (^n)  G  An,  (where 
An  is  defined  as  at  (l.f)  but  for  Ef  martingales),  such  that  the  following  three  conditions 
hold: 

(1.6)  supsupfs  [|£n(s)|]  <  oo, 

n  a<T 

(1.7)  supsup£[|^n(s)|]  <  oo, 

n  $<T 

(1.8)  lim£  [\Us)  -  f  (*«(*))!]  =  limE  [Mn (s)  -  h  (Xn(t))\)  =  0. 

n | oo  n ]oo 

Then 

(a)  Each  limit  point  of  {A„}  is  a  solution  of  the  Z?£[0,oo)  martingale  problem  for  [A,v). 

(b)  If,  in  addition,  we  assume  that  the  Z?£[0,  oo)  martingale  problem  for{A,v)  has  at  most 
one  solution,  then  Xn  =>  X ,  as  n  — ►  oo,  where  X  is  the  unique  solution  in  £)£[0,oo) 
of  the  martingale  problem  for  ( A ,  u). 

Proof  The  proof  that  each  limit  point  of  {A-,,}  is  a  solution  of  the  martingale  problem  for 
(A,  v)  will  be  analogous  to  the  proof  of  Theorem  4.8.10  [5]. 

Let  Y  be  a  limit  point  of  {A„}.  Let  (/,  h)  G  A  and  T  >  0,  and  let  (£„, <pn)  satisfy  condi¬ 
tions  (1.6  -  1.8).  Let  k  >  0,  0  <  <t2<...<tk<t<t  +  s<T.  Since  (^n,ir?n)  €  An,  it 

follows  that 

(1.9)  £[(U<  +  *)-U0-  /<+V„(u)<fu)n/l,(An(tt))|  =0, 

. v  Jt  '  .=i 
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for  all  hi,..., hk  G  B(E ),  n  £  Af,  which  implies  that 

E  [(Cn(t  +  s)-f  (Xn(t  +  s))  -  ({„(*)  -  /  (X„(t)))  -  f%n{u)  ~  h  (X„(tO)  du'j  j[h,  (Xn(t,)) 

+  E  |(/(X„(*  +  s))  -f(Xn(t))  -  j'+Sh(Xn(u))  du'j  nM*„(*0)  =  0,  Vn. 

Taking  into  account  condition  (1.8),  and  that  hi,...,hk  G  B(E),  we  have  that 

(1.10)  YimE  (. f(Xn(t  +  s))-  f(Xn(t))-  £+\(Xn(u))  du)l[hi{Xn(t ,))  =0 
and  hence 

(1.11)  E  (f(Y(l  +  s))-!(Y(t))- J'*'h(Y(u))iu)i\hAY(U))  =0. 

Thus  we  have  that  Y  is  a  solution  of  the  martingale  problem  for  (A,  v),  and  part  (a)  is 
proved. 

By  (a)  and  (b)  we  have  that  the  martingale  problem  for  ( A,v )  has  a  unique  solution. 
This  means  that  the  finite-dimensional  distributions  of  all  limit  points  of  {X„}  coincide. 
Since  the  Borel  and  cylinder  <r-algebras  are  equivalent,  the  finite-dimensional  distributions 
determine  the  probability  measure  on  I>£[0,oo).  By  this  argument,  together  with  relative 
compactness,  we  get  that  each  subsequence  of  {Xn}  contains  a  further  subsequence  which 
converges  to  the  unique  solution  of  the  martingale  problem  for  (A,i/).  By  Theorem  2.3  [2] 
we  are  done.  ■ 

In  order  to  apply  (b)  of  Theorem  1.3  we  need  to  know  how  to  determine  uniqueness  for 
solutions  of  martingale  problems  for  processes  with  values  in  arbitrary  topological  spaces. 

Theorem  1.4  Let  ( E,t )  be  an  arbitrary  topological  space,  and  let  A  C  B(E)  x  B(E). 
Suppose  that  for  each  p  €  V  {E)  any  two  solutions  X,  Y  of  the  martingale  problem  for 
(A,p)  have  the  same  one- dimensional  distributions;  i.e.  for  each  t  >  0  , 

(1.12)  P{X{t)£T}  =  P{Y{t)£T),  T  €  B(E). 

Then  any  two  solutions  of  the  martingale  problem  for(A,p)  have  the  same  finite-dimensional 
distributions;  i.e.  (1-12)  implies  uniqueness  on  the  cylinder  ty-algebra.  If  X  is  a  solution  of 
the  martingale  problem  for  (A,p)  with  respect  to  the  filtration  ,  then 

(1.13)  E  [f  (X(s  +  <))  \Tt)  =  E[f  (X(s  +  <))  | X(s)) 
for  all  f  €  B(E )  and  s,t  >  0. 

Proof  The  proof  is  completely  analogous  to  the  proof  of  Theorem  4.4.2  of  [5]. 

Corollary  1.5  Let  ( E,t )  be  a  completely  regular  topological  space  with  equivalent  Borel  and 
cylinder  a-algebras  on  £>e[0,  oo).  Let  A  G  B(E )  x  B(E).  Suppose  that  for  each  p  G  V  (E) 
any  two  solutions  X,  Y  of  the  martingale  problem  for  ( A,p )  with  sample  paths  in  De[®,  oo) 
satisfy  (1.12)  for  each  t  >  0.  Then,  for  each  p  G  V{E),  any  two  solutions  of  the  martingale 
problem  for  ( A,p )  with  sample  paths  in  Z?£[0,oo)  have  the  same  distributions  on  £>£;[0,  oo). 


4 


Proof  The  fact  that  X  and  Y  have  the  same  finite- dimensional  distributions  is  a  part  of 
Theorem  1.4.  Since  the  Borel  and  cylinder  cr-algebras  are  equivalent  the  finite-dimensional 
distributions  of  X  and  Y  determine  their  distributions  on  Z)£[0,  oo),  and  we  are  done.  | 

Let  MF(Rd)  denote  the  space  of  finite  measures  on  ( Rd ,  B(Rd))  endowed  with  the  topology 
of  weak  convergence.  For  p.  G  MF(Rd)  and  /  G  C(Rd),  let 

(114)  (f,p)  =  Jfdp. 

Denote  by  S(Rd)  the  (Schwartz)  space  of  rapidly  decreasing  functions  on  Rd,  and  by 
S'(Rd)  the  topological  dual  of  S(Rd),  the  space  of  tempered  distributions.  We  endow  S'(Rd) 
with  the  strong  topology. 

For  reasons  that  will  become  clear  later  we  now  need  to  study  stochastic  processes  taking 
values  in  MF(Rd)  x  S{Rd).  For  simplicity  denote  MF(Rd)  x  S'(Rd)  by  MF  x  S'.  Recall  that 
in  the  theory  of  weak  convergence  of  processes  in  £)£[0,oo),  where  E  is  complete  separable 
metric  space,  the  basic  equivalence  (1.5)  between  the  Borel  and  cylindrical  cr-algebras  is  of 
crucial  importance.  It  is  not,  however,  clear  that  it  carries  over  to  the  case  of  E  =  MF  x  S', 
which  is  what  we  shall  require.  The  general  problem  of  determining  which  spaces  E  do  possess 
equivalent  Borel  and  cylindrical  cr-algebras  has  been  studied  in  some  detail  by  Jakubowski  [7], 
Following  his  techniques,  we  shall  establish 

Theorem  1.6  The  space  Dmfx5'[0i  °°)  has  equivalent  Borel  and  cylinder  o -algebras. 

Before  the  proof  of  this  Theorem  we  note  the  following  lemma  from  [7]. 

Lemma  1.7  (i)  Suppose  that  the  completely  regular  topological  space  ( E,t )  has  the  fol¬ 

lowing  two  properties: 

(1)  Compact  subsets  of  E  are  metrizable. 

(2)  There  exists  a  sequence  of  {An}  of  compact  subsets  of  E  such  that  for  every 
x  €  £>£[0,1]  one  can  find  Kn  containing  the  set  x  =  {x(t)|t  G  [0, 1]}. 

Then 

(1-15)  G(DE[0,\])  =  B(DE[0,l}). 

(ii)  If  a  completely  regular  linear  topological  space  E  satisfies  (1.15)  it  also  satisfies  (1.5). 

Proof  of  Theorem  1.6  First  let  us  prove  that  MF  x  S'  generates  equivalent  Borel  and 
cylinder  cr-algebras  for  DmfXs>[0,  1], 

We  know  that  MF  is  a  separable  metric  space,  hence  it  is  homeomorphic  to  a  subset  of  R°°. 
The  space  R  has  the  properties  (1)  and  (2)  from  Lemma  1.7,  and  by  Proposition  5.3  [7]  S'  also 
satisfies  these  conditions.  Hence  by  Corollary  2.8  [7]  R°°  x  S'  has  the  property  (1.15)  and 
by  Theorem  2.1  [7]  A  x  S'  also  has  this  property  for  each  subset  A  of  R°°.  By  Jakubowski ’s 
Theorem  1.3  the  topology  on  ^£[0,1]  depends  only  on  the  topology  t  on  E,  consequently 
the  property  (1.15)  is  preserved  under  homeomorphism.  Thus,  by  a  homeomorphism  of  MF 
to  some  subset  of  R°° ,  we  have  that  (1.15)  is  satisfied  by  MF  x  S'  and  hence  by  part  (ii)  of 
Lemma  1.7  DmfXs'[Q,  oo)  has  equivalent  Borel  and  cylinder  cr-algebras,  as  required.  g 

We  need  some  additional  properties  of  the  space  MF  x  S'. 


Lemma  1.8  (i) 


B(MF  x  S')  =  B(Mf )  x  B(S'). 

(ii)  Let  N\  C  C(MF )  and  N2  C  C(S')  be  separating  for  V{MF)  and  V  (S')  respectively. 
Then 

N  =  {/1/2  :/i  6  JVi  Ul,  /2  6  U  1} 
is  separating  for  V  (MF  x  S'). 

Proof  (i)  It  is  sufficient  to  show  that 

(1.16)  B(Mf  x  S')  C  B(MF)  x  B(S'). 

Let  A  be  a  Borel  subset  of  MF  x  S'.  By  Proposition  5.3  [7]  we  have  that 
(1-17)  A=  (J  An(MFxI<n), 

nGAf 

where  Kn  are  compact  in  S'  and  metrizable.  Each  compact  metric  space  is  separable,  and 
Mf  is  a  separable  metric  space,  hence  by  [2],  p.225  we  have  that 

(1.18)  B(Mf  x  Kn)  =  B(MF)  x  B(I\n) ,  Vn. 

From  the  definition  of  the  relative  topology  in  MF  x  A'n  we  have 

(1.19)  B  (Mf  x  A'n)  =  <t(BD(Mfx  I<n),  B  er) 


=  {BH{MFx  Kn),  B  (E  B(Mf  x  S')}  , 

where  r  is  the  topology  on  MF  x  S',  and  the  last  equality  follows  by  [3],  Theorem  10.1.  By 
choice,  A  €  B(MF  x  S'),  and  so 


A  n  (MF  x  Kn)  6  B  (MF  x  Kn) .  DTIC  QT7ALITy  INSPECTED  3 


Thus,  by  (1.18), 

A  n  (MF  x  Kn)  e  B(MF)  x  B  ( Kn ) . 
By  the  same  arguments  as  in  (1.19)  we  get 

B(Kn)  =  {BnKn,  BeB(S')}, 

so  that 

B(Kn)  C  B(S'),  Vn, 

and 

A  n  (Mf  x  Kn)  €  B(Mf)  x  0(5'),  Vn. 


Acut..k>'l  .>-'or 


Mis  cra&i 
one  TAB  (7J 

!  U: . . :  f  *  r  <  : '  '"til  □ 


1  !  . 


i:  v 


Oist 


;l 


By  (1.17)  we  obtain 


>1  6  B{Mf)  x  B(S'), 


and  we  are  done. 

(ii)  The  proof  of  this  part  is  completely  analogous  to  that  of  Proposition  3.4.6  [5].  | 

Let  Ci(Rd)  denote  the  set  of  continuous  functions  with  limit  at  infinity.  In  general,  if  F 
is  a  set  of  functions  on  Rd ,  write  F+  for  {/  G  F  :  inf/**  f(x)  >  0}. 

Corollary  1.9  The  set  of  functions 

{F€  C(Mf  x  S')  :  F]uh  =  exp {(/i, Mi)  +  *  (/2,  ^2)}}  , 

where  f\  G  Ci(Rd)+,  and  /2  G  S(Rd),  is  separating  on  V(MF  x  S'). 

Proof  By  [4],  Theorem  3.2.6,  the  set  of  functions  {exp{(/,  •)},  /  G  C;(/?d)+}  is  separating 
on  V  ( MF ).  By  [6],  Theorem  3.2,  the  probability  law  of  X,  where  X  is  a  random  distribution, 
is  uniquely  determined  by  the  characteristic  functional 

C*(/)  =  £[exp{i  (/,*)}],  f€S(Rd), 

and  hence  the  set  of  functions  exp{i  {/2,  •)}  is  separating  on  V  (S').  The  result  then  follows 
by  Lemma  1.8,  in  spite  of  the  fact  that  1  0  {exp {(/,•)}>  /  G  C/(/?d)+},  since  there  exists 
{/«}  €  Ci(Rd)+,  such  that 

(1 .20)  bp-Jirn  exp  {{/,)}  =  1. 


2  Bisexual  branching  without  diffusion 

In  this  section  we  shall,  briefly,  describe  some  results  of  Kurtz  [5],  which  correspond  to  a 
generalization  of  the  bisexual  branching  system  described  in  the  Introduction,  but  for  which 
the  “particles”  perform  no  motion.  In  the  following  section  we  shall  extend  this  model  to 
the  one  of  interest  to  us,  but  it  is  worthwhile,  at  this  stage,  to  look  at  the  simpler  case. 

Kurtz  [5]  considered  a  system  made  up  of  two  types  of  particles.  Each  particle  lives  for 
an  exponentially  distributed  lifetime  with  parameter  Aj  or  A2,  depending  on  its  type.  If  a 
type  1  particle  dies  it  gives  rise  to  offspring  of  types  1  and  2  with  the  number  of  offspring 
distributed  as  (71,72).  Similarly,  if  a  type  2  particle  dies,  it  gives  rise  to  offspring  of  both 
types  distributed  as  (^1,^2)-  Assume  that  £[7<3]  <  00,  E  [0f]  <  00,  i  =  1,2.  Define 

mij  =  E  [7;] ,  m2j  =  E  [xff] ,  ;  =  1 , 2. 

Let  Zj  be  the  number  of  type  j  particles  alive  at  time  /,  and  set 
(2.1)  Z(l)  =  (2,(1),  Z,(!)),  t>  0. 
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Then  Z  is  is  a  two-type  Markov  branching  process.  We  shall  assume  that  the  process 
is  critical  and  mt:  >  0  for  all  i,j  =  1,2.  This  implies  that  there  exist  vectors  (^1,1/2)  and 
(^1 ,  ^2)  and  a  real  number  77  >  0  satisfying 


(2.2) 

/  Ai(mn  -  1) 

^lm12 

V  A2m2i 

A2(m22  -  1) 

(2.3) 

/  A](mn  -  1) 

Aimj2 

l  A2m2i 

A2(m22  -  1) 

(2.3)  (  Al(m"  ’  l)  Alfni2  )(  &  )  =  ]  . 

\  Ajm2i  ^2(77122  —  1)  /  \  £2  /  \  6  / 

Taking  i/i,i/2  >  0  and  and  £2  will  have  opposite  signs.  Kurtz  considers  a  sequence  of 
processes  {(Z[n^ ,  Z^)}  with  initial  population  sizes  (Z,  (0),  Z^O))  =  ([nzj],  [7722)),  where 
(zi,z2)  are  fixed.  Set 

(2.4)  A'W(I)  = 

(2.5)  ywW  .  + 


Then  both  and  Y(n)(<)  exp  {nrjt}  are  martingales  [1].  We  are  interested  in  the  limiting 

behaviour  of  X^n\t)  and  Y(n'(t)  as  n  —*  00.  As  might  be  expected  from  the  fact  that 
Y(n\t)  exp{nr]t}  is  a  martingale,  Y^(t)  converges  to  zero  as  n  — ♦  00,  and  so 


Z{n)(nt)  ^  6*(n)(0  ' 
n  ~  i/j6  -  V2Z1  ’ 


gfV)  f,AT<->(Q 


Thus  we  have  the  rather  surprising  result  that  the  limiting  behaviour  of  X^  gives  the 
limiting  behaviour  of  both  Z[n\nt)/n  and  for  t  >  0. 

The  limiting  behaviour  of  Y^  is  somewhat  more  delicate,  and  is  best  described  in  terms 
of  a  new  process  defined  by 

W(")m  =  Y<">(0  +  ltnr,Y^n)(s)ds. 

Jo 

To  state  Kurtz’s  main  theorem  we  require  the  random  variables 

7i  =  v\il\  ~  1)  +  "272  >  72  =  6(7i  “  1)  +  6 72, 

(2.6) 

01  =  1701  +  O2(0 2  -  1)  »  02  =  601  +  6(02  ~  1)- 

and  the  parameters  a-;  =  £’[7,7^]  and  a2;  =  £[0,-0,].  Then  Kurtz  [8]  proved 

Theorem  2.1  (a)  T7ie  sequence  {(A*”),  VY^)}  converges  in  distribution  to  a  R2-valued 
diffusion  ( X ,  W)  with  generator 

(2.7)  Af(x,  w)  =  |  (a„/IX(z,  w)  +  2aufx<w(x,  to)  +  anfww{x,  to)) 
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where  atj  =  (Ai&a-j  -  A2fia?j)/(I/i6  -  6^2). 

(b)  For  each  T  >  0, 

sup  |K^(t)  —  K^O)  exp{— nr)t}\  — >  0,  in  probability,  as  n  —*  00. 
t<T 

Consequently,  for  0  <  <1  <  t2,  it  follows  that  //*  nrjY^(s)  ds  converges  in  distribution  to 
W(t2)  —  W(<i)  as  n  — ♦  00. 

Our  aim  now  is  to  add  diffusion  into  the  above  model,  and  obtain  a  version  of  Theorem 
2.1  incorporating  the  spatial  motion. 

3  Bisexual  branching  with  diffusion 

We  start  with  a  “typical”  situation.  Assume  that  at  time  t  =  0  we  have  a  number  of  particles 
of  types  1  and  2  scattered  throughout  Rd,  which  begin  diffusing  as  independent  Brownian 
motions  with  generators  |A.  Each  particle  of  type  i,  i  =  1,2,  dies,  independently  of  the 
others,  after  an  exponential  time  with  parameter  A,.  Let  p'kl  be  the  probability  that  at  death 
a  particle  of  type  i  produces  k  offspring  of  type  1  and  l  offspring  of  type  2.  For  Borel  A  C  Rd, 
set 


Zi(t,A)  =  Number  of  particles  of  type  i  in  the  set  A  at  time  t. 


Then  Z,(t)  6  Mp  for  all  t  >  0.  Consider  the  vector  measure- valued  process 

Z(t)  =  (Z1(t),Z2(t)), 

adapted  to  the  filtration 

Tt  =  <t  (Zs  :  s  <  tj  . 

We  write 


(3.1) 


\Zi(t)\  =  Zi(t,Rd),  i  =  1,2, 


to  denote  the  total  mass  processes,  essentially  equivalent  to  the  population  processes  dis¬ 
cussed  in  the  previous  section. 

Take  now  a  sequence  Z^n\t)  =  ( Z[n\t ),  Z^lft))  of  such  processes  with  death  inten¬ 
sities  A^"\  A<n) ,  and  offspring  distributions  p\ jn\  pk\n\  adapted  to  the  filiations  Ff  = 
a  :  s  <t). 


2(") 

Pki 


Let  (71 M  ,72<">)t  have  joint  distribution  p£jn\  and  have  joint  distribution 

Set 


m 


(") 

i; 
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Let  (i/{n\  v^)T ,  (^in\^2n')T  denote  eigenvectors  corresponding  to  the  eigenvalues  r?[n) , 
—1)^  of  the  matrix  A*"',  where 


(3.2) 


A(n)  _ 


(  AfVS?  -  1)  A("W">  \ 

Af’m';’  A}l,(mS)-l) 


/ 


Then  we  make  the  following  regularity  assumptions: 


(3.3)  0  <  mj"*  <  1,  Vn,  0  <  m,y  =  lim  m)"J  <  1,  i,j  =  1,2, 

)(n) 

(3.4)  A,n  >0,  i  =  1,2,  Vn,  either  lim— f-r  or  lim-T-r-  exists,  limA-7^  =  oo,  i  =  1.2. 

n!°°  Aj  nToo  A^"J  "Too 


(") 


(3.5)  lim  —  l)  —  l)  —  m^m^]  =  0, 

A(in)A^)  ((m^  -  l)  (mffi  -  l)  - 


(3.6)  lim 

njoo 


A|n)(l  -  m^)  +  A^n,(l  -  TTI22) 

(3.7)  sup  E  [(7j(n))3]  <  00,  sup  E  [(^)3  <  00, 

A<") 

(3.8)  sup  — —  <00,  i  =  1,2. 

n  n 


<  00, 


Simple  but  tedious  calculations  show  that 


(3.9) 

-i!”> 

(3.10) 

oi"’ 

Set 

(3.11) 

A[n)A^  ((mi?  -  l)  (m&}  -  l)  -  m^m^) 

Aln)(l-m^)  +  AW(l-m^) 

A(in)(l  -  ml"’)  +  A^n)(l  -  m^)  +  0{rj[n)). 


T}\  =  limj?{n). 

n  1 00 


Under  conditions  (3.3)  -  (3.6)  we  obtain  that  |>7i|  <  00. 

Despite  the  heavy  notation  above,  things  are  are  not  as  difficult  as  they  seem.  For 
example,  if  we  return  to  the  example  of  the  Introduction,  then  we  have  A[n*  =  A =  n 
and  m,-"*  =  1/2,  i,j  =  1,2.  It  then  follows  that  j/jn*  =  0,  =  n  and  we  can  take 

j/j"*  =  i/j"*  =  1)  and  d"*  =  —  1,  £2"*  =  1  for  each  n  >  0. 

The  following  technical  lemma  follows  via  straightforward  algebra. 
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Lemma  3.1  Under  conditions  (3.3)  -  (3.6)  there  exists  a  finite  N  >  1  and  eigenvectors 
(i/1(n)  ,i 4n)),  (£,(n)  ,^n))  such  that 


(3.12) 


(3.13) 


(3.14) 


inf  u-'  >  0, 
n>N  ' 


dn)  <  0i  dn)  >  0,  v  n  >  N,  lim  |^t(n)|  <  oo, 

n Joo  I  1 


lim  L(n)£<n)  -  4n)dn)]  >  0. 

njoo 


Remark  3.2  Without  loss  of  generality  we  shall,  henceforth,  assume  that  N  —  1  in  (3.12). 

As  in  the  case  of  Kurtz’s  simpler  model,  we  need  to  define  a  few  more  additional  random 
variables  and  coefficients  before  we  can  state  our  main  result: 


(3.15) 


7"  2  -  1)  +  .  7?  =  {!”(7iw  -  1)  +  d")72,">, 

ft  s  - 1) ,  li;  s  +  d"  W  - 1). 


(3.16) 


_  W«!i”  -  Wo?!" 
~  » (d">d”  -  W) 


,*(*»)  ~ 


,  a.)  =  ,i,j  =  1,2. 


Assume  that  lima,^  <  oo  exists  (this  is  not  a  strong  restriction  given  the  results  of  the 

n|oo  1 

previous  lemma  and  assumptions  (3.3)  -  (  3.8))  and  define 


(3.17) 


We  are  finally  in  a  position  to  define  the  three  measure  valued  processes  that  interest  us. 
(Note  that  the  second  two  of  these  are  signed  measures,  or  distributions.) 


(3.18) 


n 

yW(()  _  d’,,z,(,|«  )  +  (W(0| 

WM(I)  s  yl-Ht)  +  j\^Y^\s)is  , 


In  general,  let  T>(A)  denote  the  domain  of  an  operator  A.  Here  is  th°  main  result  of  this 
paper: 
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Theorem  3.3  Let  {(X^O),  W  '"*(0))}  have  a  limiting  distribution  u  6  T\Mp  x  S),  and 
assume  that  {!^n)(0)|2}  is  uniformly  integrable.  Suppose  that  (3.3)  -  (  3.8)  hold  and 
that  limn_0oaj”)  exists  and  is  finite.  Then  (.Y(n),W(n))  =>  (X,W),  where  (X,W)  6 
CmfxS'{^,  oo)  the  unique  solution  in  DmfxS'[ 0,  oo)  of  the  following  martingale  problem 
for  (A,  i/): 

A  =  {  exp  {—  {gum)  +  i  (52,^2)}, 

(3.19)  exp  {-(gum)  +  1' (52,^2)}  (l  (-A^,  -2r?,5,  +  aug\  -  2iaugig2  -  a22g\)  ,m)  '■ 
gieT>(\A)+,g2eS(Rd)}, 

where  atj  and  g  1  are  as  in  (3.17)  and  (3.11)  and  T>(\ A)+  =  2?(|A)  n  Ci(Rd)+  . 

Furthermore,  for  each  T  >  0  and  g2  £  S(Rd), 

(3.20)  sup |^2,  ^(n)(0)  -  (52,i/(n)(0))exp{-'?2n)<}|  -+  0,  in  probability , 
as  n  —*  00. 

The  rest  of  this  paper  is  devoted  to  the  proof  of  above  theorem.  The  proof  will  rely  on 
checking  the  conditions  of  Theorem  1.3. 

4  Proofs 

There  are  a  number  of  steps  involved  in  applying  Theorem  1.3  in  order  to  prove  Theo¬ 
rem  3.3.  The  first  of  these  lies  in  finding  the  approximations  £n  and  <pn  to  f(X^n\\Vtn’>) 
and  h(X^n\W^),  where  (/,  h)  €  A.  We  do  this  in  the  following  subsection,  denoting  the 
approximations,  in  Lemma  4.1,  by  f^  and  h^K 

Section  4.2  contains  a  sequence  of  preparatory  lemmas  that  ultimately  show  that  (X^n\  W7^) 
satisfies  a  compact  containment  condition,  required  in  Section  4.3,  for  showing  that  (X^n\  lV^n^) 
is  relatively  compact  in  DmfXs'[0, 00),  a  necessary  requirement  of  Theorem  1.3. 

The  final  requirement  of  Theorem  1.3,  that  the  martingale  problem  (A,  1/)  have  a  unique 
solution,  is  established  in  Section  4.5  .  This  proof,  itself,  relies  on  a  uniqueness  result  for  a 
particular  non-linear  evolution  equation,  whose  proof  is  relegated  to  an  Appendix. 

In  Section  4.4  we  show  that  the  limit  process,  which,  according  to  Theorem  3.3  and  the 
previous  steps  is  in  DmfXs' [0, 00),  is  actually  continuous,  and  that  (3.20)  holds. 

4.1  The  martingale  approximation 

We  need  some  simple  notation  which  will  be  used  later.  If  for  sequence  of  random  variables 
{£„}  there  exists  a  constant  C  >  0,  independent  of  n,  and  another  sequence  {rpn}  such  that 

(4.1)  M<C\rPn\,  Vn 

we  shall  write  that  =  (9(0„).  Denote  by  o(l)  a  sequence  of  uniformly  bounded  random 
variables  which  converges  to  zero  as  n  — >  00. 
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Define  L  by: 


(4.2)  l  =  {/  €  C“  ( R,  x  R) :  jim  =  0,  V  t,  m  >  o}  . 

For  simplicity  we  shall  denote  by  /* ,  ,  /„ ,  fww ,  /xu,  the  first  and  second  order  partial 

derivatives  of  /  £  L. 

For  each  /  €  X,  £  C(Rd)+  (  where  C(Rd)  is  the  set  of  bounded  continuous  functions 
on  Rd  ),  and  <72  €  S(Rd)  define 

(4-3)  GLgug2  :  Mf  x  S'  — ♦  R 

by 

(4-4)  31,32(01,  02)  =  f  ((0i,0i)  ,  (02,02)),  V(/zi  ,  fj.2)  £  MF  X  S'. 

Furthermore,  define  the  operator  £  C{MF  x  S')  x  C(MF  x  S')  with  domain 

(4.5)  V(A1)  =  f£l,  gx£V(\A)+,g2eS(Rd)} 

by 

(4.6)  A\Gj,gug2{ni,ii2)  =  fz  ({0,,  0,),  (g2,  H2))  (0i0i,  01) 

“I”  2^XZ  ((01>^l)  5  (02,02))  (<*11  (0l)2 ,0l) 


+  /ru,  ((01,0l)  ,(02,02))  (a120102,  01 ) 


+  5/wu,  ((01,0l)  ,(02,02>)  (<*22  (02)2,0l) 

+  /x((0!,0 1)  ,  (02,02»  (^A0i,0i)  , 

where  the  a,j  are  defined  by  (3.17). 

Lemma  4.1  For  a//  Gs,gugi  £  V(Ai)  there  exists  (/„,  An)  €  .4„,  sucA  <Aa<  ; 

/»(<)  =  G/,„,„(X'"',|vW)  +  o  ((,<"')'’ |^M(0|), 

MO  =  +  o  (B-'|X<*|(0I)  +  o  ((,S’»)"'  |A-W(0l!) 

+0  ((,;"»)■' ixW(i)i)+ <,(1). 
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Proof  Set 


/(B)(0  =  /  ((*,  *<">(*)),  (g2,  WM(t)))  . 

To  find  so  that  (/(n>,  A(n))  €  An  ,  we  calculate  limt_0  e'1E  [/(n)(t  +  t)  -  /(n)(0 1-^?]  ant* 
obtain 

(4.7)  *W(t)  =  ( |F,(n)  ((5l,  X (52,  W(n)(0» 

+  Ain)£r  [/  ((51,  XM(t))  +  £gu(git  W^(t))  +  % g2 ) 

-  /  ((*, xwM).  (*,  W'M(O))  ] ,  4n)(0 ) 

+  (if2(n)  (<^1, AT(n)(0),  (i/2,  W/(")(0» 

+  [/  (($„*<">(<))  +  ^i,(sa,W*">(0)  +  ^2) 

-  f  ((*„  *W(«)),  ( 92 ,  W(n)(0»  ] ,  4n)(0  ) 


where 

(4.8) 


(4.9) 
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+  fww  (x,  w)  5Z(<9>52(-))2 

Jn) 

+  fx(x,  W)  -A^-) 

An) 

+  fw(x,  U>)  — ^2(-)- 

n 

Expanding  /  in  a  Taylor  series  about  (($,,  X <“>(()>,  (92. Hd"> (())),  we  have 
(4.10)  =  (If'W  «9„XI"I(()),(Sj,H'W(())) 

+  ^  [/■  «j„  X<->(()),  (ft,  H'("'(0))  £  Wl  9i 

+  /„  «J.,  *<">(<)>.  (92,  E  K]  92 

+  ~f„  «9.,A-<"'(0),(92,1V<”'(0»  £  [«?)’]  (9,)2 

+  ((9,,  XW(<)>,  (92,  »"">(<)»  «  [tt.")’]  fe)’ 

+  i/„  «Ji,Xi"l(()),(9!,lV<’»(l)»  £|7,n7;]9.92] ,  Z!”(t)) 

+  (1F2«*>«S„X<"'(()>,<9„H'<")(()» 

+  ^  [/,  «9.,  X(">(<)}.  <92.  E  [*]  9, 

+  /.  ((9.,  X  <’>(<)>.  (92,  »•*•>(«)»  £  [«]  92 
+  ((9r.X("»(i)>,(92,lVw(0))  E  [(«)’]  (91  )2 

+  i/..  ((ft,Xw(<)).(»2,H'w(0))  E  [(&)’]  (92)’ 

+  £/„  ((j.,X<">0)).  (ft.  "*•'(»)»  E  [««]  9,92]  ,  4”’(<)) 

+  /.  ((j2,X<»(<)>.  (92,  «'w(0»  (iS”  W">(f))  +  n-'O  (|X<*>(0|)  • 
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The  error  is  0(n_1|A'(n^(/)|),  since  n  l\Z\n\t)\  is  bounded  by  a  constant  times  |A"(n)(()|. 
Recalling  that 


(4.11) 


(4.12) 


7W„1  _  (0  -  ■4’-Vw(i) 

1  (  )  -  W-.W 

2  k!*^*1  -  ■'S”,dn) 


4"‘(0  = 


(*»)_(")  «r(n)  (") 

JEW!  =  -LJ-,  £W)  =  -M f- 


(4.13) 


we  see  by  (4.13)  that 


£  0»  = 


Ai"  ’ 


(4.14)  /^£K]ft,Z,W(()\  +  /^£[&]ft,4"'(<A  =  -(4">S' ■j,V'(n>(0)- 

Hence  the  terms  involving  /„,  cancel  and  (4.10)  can  be  rewritten  as 
h(n)(<)  = 

1  ((ft,-r<’)(0).(ft,W*->(0))  -rf,|f21'1  ((ft,Xw(0),(ft,IVw(0))) 

2'  kW  -  •'j-’f!'1 

+  2/.  ((ft, *<"'(<)), (ft, H^W))  lin)s. 

+  f„  ((ft , EV<”>(<)),  (ft,  W|n,(0))  «ft)  (». )2 

(4.15)  +2 /„  ((ft,A>>(0),(ft,H'M(l)))al"!lftft 

+  U  ((«„*<">(<)),  (ft,  >V<">(<)))  43*  (ft)2  ,X<">(()) 

1  /  n  (4”>F'W  ((ft, *<">(()),  (ft,  H'W(I)))  -  4”’f2,n|  «S. ,  JV«">(t)).  (j,2,  W<">(<)))) 

+  2 '  4n>{<">_4”>{|"' 

+  ((ft ,  JfW(i)),  (ft,  H""'(0))  43'  (ft  )2 

+  v~  ((«. ,  *<">(()),  (ft,  »'«■>(<)))  43'fts. 

+  /„  ((ft  ,X("’(i )),  (ft,  H,I"J(()))  43'  (ft)2 ,  >/(n)(<) )  +  0  (n-> |XW(0l)  , 


+£{ 
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where 


(4.16) 


lM  -Ain)«4n)a^  + 

'  n(W-W)  ’ 


and  a\j\  a*-n*  are  defined  in  (3.16).  Note  that,  by  3.8  and  Lemma  3.1,  the  &["*  are  uniformly 
bounded  in  n.  Finally,  by  (4.8)  and  (4.9),  we  obtain 


(4.17)i<->(<)  =  i(/,((9„XW(()>,(Sl,tVW(i)»  A*. 


+  2 /,  «a> ,  *<">«).  (p.  W<")(0))  >)!"’a. 


+ «».,  *<">(<)),  (*.  H'in,(i)»  «ft’  (s.)! 


+  2/iu-  «J„  *<">(<)).  (92,  »■'")(/)))  «<“2W 
+  /».  ((*.,  x<*>(0>.  (as,  ^'"'(O))  4'j’  te)! .  *(">(<) ) 

+  5  ( /.  ((a, ,  XW(0),  (a*  vv(">(0))  a3! 

+  /„  ((a.,  x  <■>(<)).  (*,  tv<">(0))  «i?  (a,)2 


+  2/™  ((a. .  xw«).  (9„  »•'"'(())) 

+  /„  ({a.,  x‘*>(0).<a>,»'W(t)))  4?  te)2  ,  v«">(t) )  +  o  (n-‘|A'<">(0l)  ■ 

By  (4.17)  it  is  easy  to  see  that 

(4.18)  A(n)(0  =  0  (|AT^>(0l)  - 

Set 

J(n)(i,iu)(-)  =  §(/•(*,  w)Aft(-)  +  /««(x,  to) 6i")(p1(-))2 


+  2/xu,  (x,  Itf)  fri2>^l(  )^2(  )  +  /ww  (*>  ™)  ftjSW'))  7 

and 

(4.19)  /(»>(,)  =  /W(l)  +  ((,l"))-,J<"|«s„Al">(l)),  yw(()). 

To  find  A(n)(<)  so  that  (/(n),  A(n))  €  A  we  calculate  lim«_oe_,£[/(n)(*  +  e)  -  /(nH0I^T] 
and  obtain 

(4.20)  h™(t)  =  h™(t)  +  o  ((^n))_1  l*(n)(0l2)  +  O  ((4n))~'  !*(n)(0l) 

-  (J<*>  (jJtww(t)))  ,y(n)(0)  • 
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Combining  terms  we  obtain  the  desired  result.  I 

Conditions  (1.6) —  (1.8)  now  follow  immediately  from  the  following  bound. 

Lemma  4.2  £[sup(<T  |X^n^(t)|2]  is  uniformly  bounded  in  n  for  each  T  >  0. 

Proof  From  (4.7)  we  have  that,  for  all  e  >  0, 

(4.21)  M,»>(l)s  exp  {-e|x<">(l)f} 

-  £  <%•  [exp  j-e  (|X<">  (S)|  +  £)’}  -  exp  {-<  [X<“>  (>)f }]  , 

+  [expj-<  ^X1”1  (s)|  4  |  -  exp  (-£  |X,n|  (-s)j2  j 

is  a  martingale.  Hence 
E[t|x<n»(t)|2exp{-e|A-<n>  (t)f}] 

<  £  [l  -  exp  IjV^)  (0|2}] 

=  1  -  £  [exp  (0)|2}] 

—  E  y*  (x[n)Er  [exp  {-<  (s)|  +  ^  J  -  exp  {-e  (*<»>  (s)|2}]  ,  £,(n)(^ 

+  ^n)£^|exp|-f^n)(s)|  +  ^  j-exp{-£|xH(s)|2}  ,Z[n)(s)j  ds 
=  l-£[exP{-c|^(0)|2}] 

+  £  [^(Ai">exp{-e  |*(n)w|2} 

x  Er  Ji_exp|-e^|x(")(a)|  +  ^)2-|x(n)(s)|2j|  ,z[n\s)j 
+  ^n)exp{-c|A-(n)(s)|2} 

x£*„  1-expj-e  ^|*<">(s)|+^  -|^(s)|2jj  ,Z{2n\s)^ds 

<  l-£[exp{-f|jr(n)(0)|2}] 
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+  E 

\lo  ^Ain)exp{-f  |^(n)  (*)f  }€£7n  [ 

(|xW(5)|  +  ^)2-|xW(S)2 

,z.(T°(0^ 

+  ^A5n)exp|-c|x{n)(s)|2|e^n 

^|x(n)(o|  +  ^  -|*(n)(of 

,zin)(s)\  ds 

<  tE  |X(n)(0)|2 


jr'exp{-<|xWWf}^£r[W)IJ 


z[n\s) 


n 


+  2-L-Er  [7r]|x^(S)||xi(n,(5)| 


\(n)  \ 

+  2^-^n[^]|xW(S)||4n)(5)|j  ds  . 

By  the  uniform  boundedness  of  A-n*/n,  f^n[(5")2]*  E^n [(^" )2],  and  the  definition  of  X(n>,  there 
exists  a  constant  C\ ,  independent  of  n,  such  that 

(4.22)  E  [c  |x<">  (of  exp  {-c  |X<">  (Of  }] 

<  <E  [|*<">  (0)f]  +  cC,  jf  £  [|x<n)  W|]  ds 


+  2ai[n)j*E  [exp  {-e  |*<">  (s)f}  |*<">  (s)f 


ds. 


Recalling  that  |Xf"^(t)|exp  is  a  martingale  (c.f.  [1]),  (4.22)  can  be  rewritten  as: 

(4.23)  E  [c  |X<">  (Of  «P  {  -<  |*(n)  (of  }] 

<  c£[|xW(0)|2]+€Cij[‘  exp  {^n)a}£;[|xW(0)|]  ds 
+  2aj[n^  J*E  [exp  {-c  |x<n>  (s)f}  |*<n>  (s)f]  ds 

<  tE  [|X<")(0)|2]  +  fCitexp  {»fc(n)t}£  [|*(n)  (0)|] 

+  2 '*&n)J*E  [exp  |-c  |x<n>  (s)f  }  \x™  (a)f  ]  ds. 


By  (4.23)  and  Gronwall’s  inequality  we  obtain 
(4.24)  E  [f  |x<">  (Of  exp  {  -c  |x<">  (l)f  }] 

<  €  (e  [|x(n)  (0)fj  +  C, t exp  }£  [|x(n)  (0)|] )  exp  {2^n)f} 
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Dividing  by  t  and  letting  e  1  0  gives 

(4.25)  £[|*(n)(0f]  <  [e  [|^(n)(0)|2]  +C1texp{7?1(n)<}£[|x(")(0)|j)exp{2r?jn)<}. 

That  is,  £[|*(n,(0l2]  is  uniformly  bounded  in  n .  Once  again  using  the  martingale  properties  of 
|A^n)(<)|,  we  obtain,  by  Doob’s  inequality,  that  £[supKT  |X^(£)|2]  is  also  uniformly  bounded  in 
n  for  each  T  >  0.  ■ 


4.2  Compact  containment 

Definition  4.3  Let  ( E,t )  be  a  completely  regular  topological  space  with  metrizable  compacts.  Sup¬ 
pose  that  n  =  1,2,...,  is  a  Tf -adapted  process  with  sample  paths  in  i?£[0,oo).  We  say  that 

the  compact  containment  condition  holds  for  {X^},  if,  for  every  c  >  0  and  T  >  0,  there  exists  a 
compact  set  Tt<r  C  E  for  which 

inf  P{x(">(t)  €  re,r  for  all  0  <  t  <  t}  >  1  -  c. 

In  the  following  lemmas  we  shall  denote  £ f=1  (d,g)2  by  {dg)2. 

Lemma  4.4  E[(g2,  IV(n)(<))2]  is  bounded  uniformly  in  n  for  each  g2  €  S{Rd). 

Proof  It  is  easy  to  see  that  for  each  c  >  0 

exp{-e(|^i|2  +  ($2,/i2)2)}  €  Z>(4i). 


From  (4.7)  we  obtain  that 


A//2)(t)  =  exp  | 

-c(|x(n)(i)|2  +  (g2,lFW(t))2)j 

(4.26) 

.|-€(|xW(»)|2  +  (j2,WW(a))2)} 

X 

(2c(ff2,W<»>(s))2  -  l)  (£pj  {dg2)2  -  {g2,WW{s))^&g2 

+  *[n)Er  |exp|-e^(|xW(s)|  +  ^)2+  ((<?2,^n)(s))  +  } 

-  exp  {-€  (|X(")(S)|2  +  (ff2,^")(0>2)}]  ,^n)(^)) 

+  ^cexp  {-c(|xW(s)|2  +  (ff2,iy(")(s))2)} 


(2c(ff2,VFW(s))2  -  l)  (dg2)2  -  (92,WW(s))^Ag2 


20 


+  A <%„  exp  |-f  ^|*(n)(s)|  +  &)’  +  ((ft,WW(.))  +  f  52)2j  | 

-  exp  |-f  (|*<">  (5)|2  +  (ff2,^(")(0)2)  }]  ,  z‘n)(s)} 

-2€(S2,WW(s))exp{-€  (|x<">(s)f  + 

is  a  martingale.  By  simple  calculations  (as  in  the  previous  lemma)  we  obtain  that 

E  [f(^2,^n)(t))2exp  j-e  ^|x(n>(<)|2  +  (s2,W(n)(f))2j  j] 

<  <E  [o  (|*(n)  (0)|2)]  +  ifoE  [o  (|*<n>  (s)f)]  ds  +  (J*E  [o  (|*<">  (*)|)]  ds 
+  \J*E  [«xp  j-c  (|x(n)(s)|2  +  (52,^")(S))2)|(ff2,W'(")(s)}2]  ds 

=  (Bn(t)  +  \JoE  [eexp  |-c  (|X(n>  (s)|2  +  (52,  W(n)(s)}2)  }($2,  W^H*))2  ds, 

where  the  Bn(t)  >  0  are  bounded  uniformly  in  n.  (This  follows  from  the  uniform  boundedness 
of  L’fsup,^  |Af(")(s)j2]).  Hence  there  exists  a  function  B  such  that  B(t)  >  Bn(t)  for  each  n.  By 
Gronwall’s  inequality  we  obtain 

E  [(</2,^<n)(0)2exp  {-c  (Of  +  (<72,  W(n)(0}2)}]  <  5n(0exp{l<},  Vc  >  0. 

Letting  (  |  0  we  find,  by  monotone  convergence, 

(4-27)  E  [(<72,W(n>(<))2]  <  ^n(<)exp{i<}  <  B(t)exp{i<}, 

so  that  £[(521  W^"H0)2]  is  uniformly  bounded  in  n. 

■ 

Lemma  4.5 

(4.28)  (<72,  ^(f))  -  \j‘{Ag2,Y^(s))  ds 
is  a  martingale  for  each  n  and  g2  €  S(Rd). 

Proof  Note  that  exp{-e|/iij}sin(£(02,/*2))  €  V(Ai).  Define  a  martingale  Mi3\t )  as  in  the 
previous  cases  : 

M<3>(t)  =  exp  (-£  |x<">  (t)|)  sin  (((ft,WW(l))j 

(4.29) 

~lo  ^€eXp{-f  |^(n)(S)|}  ^cos(f(?2,^n)(s)))^-Ap2 
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-tsin  (e(g2,WM(s)))  (^2)2j  ,Z<n)(0^ 

+  ^eexp{-e|A:(n)  (s)j}  ^cos  (c($2,  W(n)(5)^ ~Ag2 

-esm  (c(g2,WW(s)))  (^j  (%2)2j  ,Z<n)(^ 

+  (A [n)Er  [exp  {-e  ([*<">  (s) |  +  }  sin  (e(y2,  WW(s))  +  ^52) 

-  exp  {-£  |AT(n)  (s)|}  sin  (c(s2,  ]  ,  z{n)(s)) 

+  (\X2n)Ej«  exp|-£  (|*«(,)|  +  ^}sin  ^f(<?2,  H/(n,(s))  +  ^-g2) 

-  exp  {-e  |x<">  (s)|}  sin  (e(g2,  W  W(s)))]  ,  Z<n)(s)) 

+  cexp  {-e  |*<»>  (s)|}  cos  (<(g2,  WM(s)))r1^(g2,  Y^(s))  ds 

Let 

(4.30)  A/W(t)  =  lim 

=  (92,WW(t))  -  J*  ^Er  [£] 

+  ^4%-  [^]  92,Z<n\s)j  +  i(Aj2,yW(S)) 

+  f?2n)(ff2,y(n)(s))ds 

=  (g2,w^)(t))  -  \J*(a92,yW(s))  ds , 

where  the  last  equality  follows  by  definition  of  yW.  The  above  limit  exists  almost  surely.  On  the 
other  hand,  it  is  easy  to  see  from  (4.29)  that 

<  |(*,W*t>(0)|  +  j[<Ci|A'H(s)||(S2,iyW(s))|  +  C2|x(")(s)|  ds 

where  C\  >  0  and  C2  >  0  do  not  depend  on  c.  From  the  uniform  boundedness  of  £'[|A'(n)(t)|2]  and 
,  we  obtain  the  uniform  integrability  of  \Mi3\t)/t\,  which  means  that  convergence 
is  in  L\.  Consequently  M^\t)  is  a  martingale.  | 

Lemma  4.0  £[sup(<7-  (g2,  W(n)(t))2]  is  bounded  uniformly  in  n  for  each  g2  €  S(Rd). 
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Proof  From  the  previous  lemma  we  have  that 


(S2,w(n)(0)  =  M^(t)  +  ^^{Ag2,Y^(s))ds. 

Then 

E  sup(<72,VF(n>(0)2  =  E  |sup^A/(3)(<)+5^<(A52,K(n)(s))^) 

<  2E  jsup  (M(3)(t))2j  +  2E  sup*  (JQ  (^S2-y(n>(5))  ds'j 

<  8E  [(a/(3)(T))2]  +  \\\Ag2\\2T2CE  |sup|x(n>(0|2 

<  8 E  2(?2,"/{n)(i')}2  +  i  (a92, Y{n\s))  dsj 
+  i||A<72||2r2C£  sup|x(")(o|2  , 

where  C  is  the  constant  such  that 

|(Aff2,y(n)(0)|2  <  C||A^||2|x(")(0|2,  Vn. 


Finally  we  obtain  that 


E  sup^ff2i  S  1  §E  (g2,WW(T))  +  || Ac/2 ||2  T2CE  sup  |AT(n^  (t)| 


FYom  Lemmas  4.2  and  4.4  we  obtain  uniform  in  n  boundedness  of  £[supt<7  (g2,  ]. 


■ 


Lemma  4.7  For  each  gx  €  2>(£A)+,  g2  €  S{Rd),  the  processes  (gi ,X<n)(0)>  <ff2,  W(n)(0}  and 
((5ii-^(n,(0)>  (ff2,W(n,(0))  satisfy  the  compact  containment  condition. 

Proof  The  proof  is  immediate  by  the  uniform  boundedness  of 

E  sup(<72i  i  E  sup(^i,XW(0) 

t<T  '  '  L*<7X 

and  Chebyshev’s  inequality  in  the  case  of  the  processes  (gi,X^n\t)),  {92,W(n\t)).  In  the  case  of 
the  process  ({$,,  *<">(*)),  (g2,  W(n)(0))  the  result  follows  from  the  fact  that  the  product  of  two 
compact  sets  is  compact  in  the  product  space.  I 
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4.3  Relative  compactness 

We  now  check  the  relative  compactness  of  {(A(n), W(n))}.  En  passant,  we  shall  show  that  X 
converges  weakly  to  super-Brownian  motion;  i.e.  the  process  with  values  in  Mp  which  is  the  unique 
solution  of  the  martingale  problem  for  A*  with 

A*  =  | exp  {—  (g\,p)},  exp  {-  (gi,p)}  (~V\9i  +  ^ix  (Si)2  ~  ■  5i  €  V  (|A)  + j  • 

Recall  that  V(\A)+  =  V{\A)  n  C,(Rd)+. 

Lemma  4.8  The  following  sequences  of  processes  are  tight  for  each  g\  €  7?(^A)+,  g 2  €  S(Rd): 

1)  {(($i,*<">(0),(ff2,  WW(t)))}  m  Dr,[ 0,oo), 

2)  {(p„ *(")(<))},  {{52,WW(0)}  in0*[O,oo), 

3)  {(51,XWW)  +  (52,WW(t))}  in  Dr[ 0,oo). 


Proof  2  and  3  are  simple  corollaries  of  1,  which  we  now  prove.  By  Lemma  4.1  we  obtain  that, 
for  all  Gjt3ug2  €  V(A\),  we  can  choose  the  processes  f^  and  h^n\  (/*n\  h^)  6  An  such  that 


lim  E 

njoo 


sup  [/<”>(<)  -  /  ((<7X ,  XW(t)),  (g2,  W(")(0))  | 


lim  E 

ntoo 


supo((f£,))'Vn)(0l) 


=  0. 


The  last  line  follows  from  the  uniform  boundedness  of  £[supt<T  |A^n^(t)|].  Furthermore, 


supi? 


sup 

t<T 


|/x(n,(0| 


<  sup£ 

n 

+  sup£ 
n 

+  sup£ 


+  supE 

n 


sup  \Aj((9i,X^{t)),  (g2,W^(t)))\ 
sup  0  (n-1  |a»  (t) |) 
sup0((^n))"1|A-(")(t)|) 
supo((r?<n))“1|AH(0|2) 


+  sup  0(1) 


<  00, 

where  the  final  line  follows  from  the  uniform  boundedness  of  /?[supt<]>  |A^n)(t)|2],  and  the  bound¬ 
edness  of  functions  in  the  range  of  the  operator  At  .  Since  this  holds  for  each  choice  of  g\,  g2,  and 
by  Lemma  4.7  we  have  compact  containment,  applying  Theorems  3.9.1,  3.9.4  [5]  now  completes  the 
proof.  I 

Lemma  4.8  and  Mitoma’s  theorem  [10]  immediately  yield 

Lemma  4.9  {W(")}  is  a  tight  sequence  of  processes  in  Z?s<[0,oo). 

Lemma  4.10  {A’*")}  converges  weakly  to  super-Brownian  motion. 
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Proof  We  can  take  g2  =  0  and  thus  obtain  A  C  C(Mp)  x  C(Mp),  where  A  C  A\  is  defined  as 
follows: 

D(A)  =  linear  span  of  jexp  {- (gi,p)},  9\  €  V  +  | 

A  exp  {-  (gi,n)}  =  exp  {- (gup)}  (~Vigi  +  (5i)2  -  ■ 

By  Theorem  2.4  [14]  and  results  of  Perkins  [12]  the  closure  of  A  generates  a  strongly  continuous 
contraction  semigroup  on  C(Mp)  (with  corresponding  super- Brownian  motion  X)  and  D(A)  is  a 
core  for  A.  By  Theorem  3.2.6  [4]  the  set  of  functions  {exp  {-  strongly  separates  points  in 

Mp .  Thus,  all  the  assumptions  of  Theorem  4.8.2  [5]  are  satisfied  and  we  are  done.  | 

Lemma  4.11  {(X(n),W(n))}  is  a  relatively  compact  sequence  of  processes. 

Proof  We  check  the  conditions  of  Theorem  4.6  [7]  (Jakubowski’s  tightness  criterion).  {A'(n>} 
and  {W(n)}  are  tight.  This  implies  that  the  compact  containment  condition  holds  for  each  of  these 
processes  and  by  the  same  arguments  as  in  Lemma  4.7  it  holds  for  the  pair  {(.Ar(n),  Define 

the  family  of  functions  F  :  Mp  x  S'  — >  R  by 

=  (guni)  +  (g2iHi) ,  gi  e  02  e  S(Rd). 

This  family  separates  points  in  Mp  x  S'  and  is  closed  under  addition.  By  Lemma  4.8  the  sequence 
{/(A*"),  W^)}  is  tight  for  all  /  6  F.  Jakubowski’s  conditions  are  satisfied  and  the  proof  is 
complete.  ■ 

Theorem  1.3  and  Lemmas  4.11,  4.1  now  imply 

Corollary  4.12  Each  limit  point  of  {(A^n),  W(n))}  is  a  solution  of  the  £>A/rX5'[0, 00)  martingale 
problem  for  (Ai,v). 

4.4  Continuity  of  the  limit  process 

Throughout  this  subsection  we  shall  denote  by  W  one  of  the  limit  points  of  the  sequence 
and  all  results  are  obtained  for  each  limit  point  of  [H,W],  For  simplicity  of  notation  we  shall 
denote  (</,W(t))  and  ( g,X(t ))  by  Wt  (g)  and  Xt{g)  respectively,  while  the  increasing  process  of 
each  martingale  Mt  will  be  ( M)t . 

We  shall  treat  only  the  continuity  of  the  process  Wt ,  since  the  continuity  of  X  follows  from  the 
fact  that  X  is  super  Brownian  motion.  Using  this  result  we  shall  prove  the  last  part  of  Theorem  3.3, 
viz.  that  -  y^nH0)exP  converges  in  probability  to  zero  for  each  #2  €  S{Rd). 

Lemma  4.13  For  each  g2  €  S(Rd),  Wt  (g2)  is  a  square  integrable  martingale  with  increasing 
process 

(4.31)  (W{gi))t  =  J*a22X.(gl)  ds. 


25 


Proof  The  uniform  boundedness  in  n  of  l?[(W/n^(02))2]  (obtained  during  the  proof  of  the  compact 
containment  condition)  immediately  gives  the  existence  of  second  moments  of  Wt  (<72)-  We  now  use 
the  same  “trick”  as  in  Lemma  4.5. 

(4.32)  M^\t)  =  exp  (-( |A,|)  sin  (eWt(g2)) 

-  J  «2  exp  {-e  \X.\)  ^sin  (tW,  (92))^  \X,\ 

-  cos  {eW,  (g2))ai2Xt  (92)  ~  sin  ((W,  ($2))^y  (s2)  j 

-  e exp  {-e  IXjI}  sin  {eW a  (p2))f?i  |A,|  ds 

is  a  martingale.  Since 

(4.33)  lim  (  -K-Z  =  Wt(g2 )  a.s., 

the  dominated  convergence  theorem  gives  that  T  t  (92)  is  a  martingale.  Now  note  that 

(4.34)  M«\t)  =  exp{-£|Xf|)sin(€Wap2)2} 

-  Jo  f  exp  {-e  (f  sin  (clU,  (p2)2)  ~  |Jf,| 

-  2e  cos  (cW.  (g2)2)  W.  (g2)  auX,  (g2)  +  cos  (cW,  (52)2)a22X5  (g2) 

-  2  tW,  (g2)2  sin  (e W ,  ( g2f)a22X,  (p2) 

-  sin  ((Ws  (g2)2)v  1  E X  g | )  ds 

is  a  martingale.  Under  the  same  arguments  as  in  (4.33) 

(2) 

(4.35)  Wt(g2)2-  (  a22X,  (g2)  ds  =  limM<  — ,  a.s. 
is  a  martingale  and  so  (W  (g2))t  =  /J  a22X ,  (g2)  ds. 

» 

Lemma  4.14  Assume  that  W0  €  S'.  Then,  for  each  g2  €  S(Rd),  Wt(g2)  is  a  continuous  martin- 
gate. 

Proof  Define 

(4.36)  WtK  (ga)  =  Wt  (g2)  -  K,  for  all  K  €  R. 
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Once  again  we  use  the  same  “trick”. 

M^(t)  =  exp  {-€  M  -  (W?  ( g2 )2  -  W*  {g2f)  } 

-  J\  exp  {-( |*|  -  (W?  fa?  -  Wf  ( g2f )  }  (c2^1 1*|  +  2  eW?  (g2)2  a„*  (ff2) 

+  (2l*A  (02 )2  -  l)  022*  (ff2)  ~  l*l) 
is  a  martingale.  Under  the  same  arguments  as  in  (4.33)  we  obtain  that 
A/<3>(<)  =  exp  {-(wtK(g2)2-WQK(g2)2)} 

-  f  exp  {-  (W*  (g2)2  -  Wq  (02)2)}  (2 W*  ( g2 )2  -  l)  o„*  {g\)  ds 
JO 

—  lim a.s. 

is  a  martingale.  By  Corollary  2.3.3  from  [5]  we  have  that  the  process 
M(4>(<)  =  exp  {-  (wtK  ( 02 )2  -  W*  (02)2)  -  J‘  (2 W,K  (02)2  -  l)  o22*  (ff22)  ds } 

=  exp  {-  (lUA'  (02)2  -  W0K  (02)2)  +  jf  a 22*  (s22)  ds  -  ($2)2  *  (s22)  *} 

is,  at  least,  local  martingale.  By  setting  g2  =  m2g2  in  the  above,  we  see  that,  for  every  m  €  R, 
(4.37)  exp  j-m2  |^1UA'  (02)2  -  (02)2)  -  a22X,  (g2)  dsj 

-^/o'4a22VUA  (02)2*  (ff2)  ds} 
is  a  local  martingale.  By  Lemma  4.13 

WtK  (02 )2  -  W*  (02 )2  -  Jl a 22*  (52)  ds 

is  a  martingale.  The  continuity  of  this  martingale  derives  from  the  result  given  in  the  following 
lemma  [13]: 

Lemma  4.15  Let  X  be  a  local  martingale  such  that  Xq  =  0  and  A  be  a  continuous  increasing 
process  such  that  Ao  =  0.  //exp  {kXt  -  (k2/2)At}  is  a  local  martingale  for  every  k  6  R+  ,  then  A 
is  continuous  and  At  =  ( X)t . 

Thus,  by  (4.37),  we  have  that  WtK  (g2)2  -  Wff  (02 )2  is  a  continuous  submartingale.  This  is  true 
for  every  K  €  R  and,  thus, 

WtK(g2)2-W?(g2)2  =  (Wt(g2)-K)2-Wt(g2)2 
=  K2  -  2KWt  (02) 

is  a  continuous  process  for  each  K,  which  implies  that  Wt  (g2)  is  continuous.  | 

Lemma  4.14  and  Theorem  1  [9]  immediately  yield 
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Corollary  4.16  Wt  is  a  continuous  process  taking  values  in  S'. 


Corollary  4.17  ForT  >  0,  sup(<j- |(<72>  F^(t))  —  (<72)F(n)(0))  exP  {— converges  to  zero  in 
probability  for  each  g2  £  S(Rd). 

Proof  The  proof  is  completely  analogous  to  that  of  Theorem  9.2.1(b)  [5],  using  the  fact  that,  for 
each  02  €  S(Rd),  Wt  ( g 2)  is  a  continuous  process.  g 


4.5  Uniqueness 

In  this  subsection  we  shall  prove  that  the  martingale  problem  for  A,  introduced  in  Section  3  as: 

(4.38)  A  =  {  exp  {-  (guHi)  +  i  (92,1*2)}, 

exp{-  (gi, Mi)  +  i(g2,P2)}(\  (-Affi  ~  2r?iSi  +  anff?  -  2ia12gxg2  -  a22g2)  ,pi)  ■ 
9leV(\A)+,  g2eS(Rd)}, 

has  a  unique  solution. 

Throughout  this  article  we  worked  with  real-valued  functions.  Now  for  simplicity,  we  switch  to 
complex- valued  functions.  We  obtained  in  Corollary  4.12  that  all  the  limit  points  of  {(X^n\H^n))} 
are  solutions  of  the  martingale  problem  for  (A\,  u).  By  the  obvious  fact  that 

{  exP{-(?1,^1)}sin((ff2,/X2))Uexp{-(ffi,/i1)}cos((p2,^2»  ■  9\  €  V{A)+  ,g2  £  s}  C  2>(Xi), 

we  obtain  that  all  the  limit  points  are  also  solutions  for  ( A ,  1/).  Thus,  in  order  to  prove  weak 
convergence,  it  is  sufficient  to  prove  uniqueness  for  ( A ,  u).  Let  (X,  W)  be  any  limit  point  of 
W(n))}.  For  simplicity,  in  the  next  two  lemmas  we  shall  use  the  notation  introduced  in 
Section  4.4;  that  is  (g,W(t))  and  ( g,X(t ))  will  be  denoted  by  Wt(g)  and  Xt(g),  respectively. 

Lemma  4.18  For  each  g\  £  Z>(jA)+  ,  g2  £  S{Rd),  the  bracket  process  of  X(  (g  1)  and  Wt  ( g2 )  is 

(4.39)  WffO.W'fo)),  =  [tXa(aug1g2)ds. 

Jo 

Proof  By  Ito’s  formula 

exp{-ATj  (ffi)  +  iWt  (g2)}  = 

exp{-X0(g\)  + iW0(g2)}  +  exp  {-A,  ($i)  +  iW,  (g2)}Xs  (-pigi  -  5A01)  ds 

+  \  /‘exp  {-X,  (ffl)  +  iW.  (*)}<*«-*(*,)  +  iW(g2))s)  +  Mu 
Jo 

where  M*  is,  at  least,  a  local  martingale.  Recall  that 
expf-A',  (gi)  +  iWt(g2)} 

-  ijf  exp  {-A,  ( gi )  +  iW,  (g2)}Xt  (-2^0!  -  Affi  +  an(3i)2  -  2ial2gig2  -  a22(g2)2^j  ds 


28 


is  a  martingale  and 


(X(g))i  =  £x3(ang2)ds,  V5€P(  §A), 

(W(g))t  =  £x,  (a22g2)  ds,  Vg  £  S(Rd). 

Combining  the  above  result  with  the  fact  that 

{-X(gi)  +  iW(g2))t  =  (X  (gi))t  -  (W(g2))t  -  2t  {X  (gi),W(g2)),  , 

we  are  done.  ■ 

Let  Co(Rd)  be  the  set  of  continuous  functions  tending  to  zero  at  infinity.  Introduce  the  additional 
notation 

(4.40)  C'(Rd)  =  {g:  g  =  gi  +  ig2-.gieCl(Rd),  g2eC0(Rd)}, 

(4.41)  C'(Rd)+  =  {g:  geC'(Rd),  gi€V(±A)+,  g2ev('iA)nCo(Rd)y 
and  the  operator 

(4.42)  A  -  {  exp  {-(guPi)  +  i  (92^2)}, 

exp{-  (51, Hi)  +  i(g2,>i 2)}  (5  (-Ajf!  -2Tjigi  +ang2  -2iaugxg2  -  a22g])  ,m)  : 

Then  we  have 

Lemma  4.19  ( X ,  W)  is  a  solution  of  the  martingale  problem  for  (4,  v)  if,  and  only  if,  it  is  a 
solution  of  the  martingale  problem  for  (A,  u). 

Proof  Since  A  C  A,  if  (X,  W)  is  a  solution  for  (A,  v),  then  it  is  also  a  solution  for  (4,  v ).  Let 
(A",  W)  be  a  solution  for  (A,  ts).  For  each  gi  £  C\Rd)+  ,  g2  €  S{Rd),  we  obtain,  by  Ito’s  formula, 
that 


exp{-Xt(g\)  +  iWt(g2)} 

=  exp{-X0(gi)  +  iW0(g2)}  +  exp  {-X,  (gi)  +  iWt  (g2)}Xt  -  ^Affi)  ds 

+  i  f\xp{-X.(9l)  +  iWt(g2)}d((-X(9l)  +  iW(g2))t)  +  Mu 
Jo 

=  exp  {-Xo(gi)  +  iWo(g2))  +  \  exp  {-A,  ($,)  +  (y2)} 

x  X,  {-2t)igi  -  Aji  +  an(ffi)2  -  2iaugi92  -  022(02)*)  ds  +  Mt , 
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where  the  second  equality  follows  by  previous  lemma,  and  Mt  is,  at  least,  a  local  martingale.  But 
all  the  terms  in  the  right  and  left  hand  sides  are  bounded,  so  that  Mt  is  a  martingale.  g 

In  the  proof  of  uniqueness  the  following  lemma  will  play  a  crucial  role: 


Lemma  4.20  The  nonlinear  evolution  equation 


(4.43) 


dU(t) 

dt 

.  U{  0) 


£/(t)2  +  U(t)rfi  +  ^A U{t)  +  ia\2U{t)g2  +  \o.22  (92 )2 


9\ , 


where  gi  £  C'(Rd)+  ,  g2  €  S(Rd),  has  a  unique  strong  solution  on  R+  ,  such  that  U(t),  dU  (t)  /dt, 
and  5 A U(t)  are  continuous  functions  from  R+  to  C'(Rd),  and  U(t)  £  C'(Rd)+  for  each  t  £  R+  . 


Proof  Appendix. 

Lemma  4.21  If  (X(t),W(t))  is  a  solution  of  the  martingale  problem  for  (A,  u)  in  Z?AffxS'[0,  oc), 
then,  for  each  gx  £  C'(Rd)+,  T  >  0, 

(4-44)  exp  {-  (U(T  ~  t ),  X(t))  +  i  (g2,  W(0>} 

is  a  martingale  for  0  <  t  <  T  ,  where  U(t)  is  given  by  the  unique  solution  of  (4-43)  with 

u(o)  =  giec\Rd)+. 


Proof  Define 

(4-45)  F(9i,g2)  =  Vi9i  -  jflxi  (ffi)2  +  ian9i92  +  §<*22  U2)2  • 

Let  X(t)  =  (.Y(t),  W(t))  be  a  solution  of  the  martingale  problem  for  (A,  1/),  v  £  V(Mp  x  S'). 
Define 

(4.46)  «(*,*(*))  =  exp{-(U(T-s),X(t))  +  i(g2,W(t))},  VO  <  s  <  <  <  7\ 

For  each  (#*i, p2)  €  Mp  x  S'  we  have 

du(s,(pup2)) 


(4.47) 

Hence 


ds 


=  exp  {-  (U(T  -  s), /zi)  +  t(ff2,M2>} 


dU  (s) 


ds 


i=T-s 


-Ml 


E  (t2,X(t2fj  -  u  (fi,AT(t2))| 


=  E 


It 


”  exp  {-  (U(T  -  s),X(t2))  +  t  (g2,W(t2))}  ^  dU  {S)- 


ds 


8 


,X(t2))ds 


=  Ey\(s,X(t2j)  ds  *,]. 


Fix  0  <  t\  <T.  Since  g\  £  C\Rd)+,  we  have  U(T-t\)  £  C\Rd)+,  for  every  0  <  ti  <  T.  Therefore, 
by  the  definition  of  the  martingale  problem  for  (A,  u),  we  have  that 

exp  {-  ( U(T  -  h)9i,X(t))  +  i(g2,W(t))} 

-  J\xp{-  (U(T  -  h),X(s))  +  *(52,^))}  x  {-F(U(T  -  U),g2)  -  lA U(T  -  U),X(s))  ds 
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is  a  martingale  for  0  <  t  <  T.  Define  the  function  tr  by  the  following: 
(4.48)  £;[«(t1,X(t2))-«(t1,X(t1))|/-1] 

=  E  [exp  {-  (U(T  -  ti),  X{t2))  +  i  (g2,  W{i2))} 


-  exp  {-  (U(T  -  M,X(<i)>  +  i  (p2,  W(t,))}|  T tj 
=  E  [I*  exp  {-  (U(T  -  ty ), X(s))  +  i  (g2,  W(s))} 

x  (-F(U(T~tr),g2)-  \AV{T-U),X(sj)  ds|jF(l] 


=  E  w  ^i,s,  X(s))  ds 


,  VO  <  tj  <  <2  <  T  <oo 


It  is  easy  to  check  that  all  conditions  of  the  Theorem  4.3.4  [5]  are  satisfied  and  hence 


(4.49) 


u(«,X(<))  -  {t>(s,X(s))  +w  (s,s,X(s))}  ds 


is  an  Tt -martingale.  Finally,  substituting  the  definitions  of  u,  v,  w  we  obtain  the  desired  result.  § 


Theorem  4.22  The  martingale  problem  for  (A,  v)  has  a  unique  solution. 

Proof  Let  ( X(t),W(t ))  be  any  solution  of  the  martingale  problem  for  (A,i/),  v  €  V(Mf  x  5'). 
Take  g2  €  S(Rd)  and  gi  €  V(^A)+  .  (Recall  this  means  that  gi  =  gn  +  igu  €  C\Rd)+  ,  gn  =  0.) 
Then  by  the  previous  lemma,  setting  T  =  t,  we  obtain 

£[exp{-  (gx  ,X(t))  +  i(g2,lF(t))}]  =  £[exp{-  (U(t),X( 0))  +  i(g2,W(0))}],  VO  <  t  <  oo, 

where  U(t)  is  the  unique  solution  of  (4.43)  with  17 (0)  =  g\.  Thus,  by  Corollary  1.9,  we  obtain 
that  that  any  two  solutions  of  the  martingale  problem  for  (A,  v)  have  the  same  one-dimensional 
distributions.  By  Theorem  1.4  the  desired  result  follows.  ■ 

Remark  4.23  The  proof  of  Theorem  3.3  is  now  finished,  since,  by  Corollary  4.12  and  Theo¬ 
rem  4-32,  all  the  conditions  of  Theorem  1.3  hold. 


Appendix:  The  solution  of  a  particular  non-linear  evo 
lution  equation 


Lemma  A.l  The  nonlinear  evolution  equation 


dU[t) 


m 


anU(t)2  +  U(t)rjr  +  \AU(t)  +  ial2U{t)g2  +  Ia22  (g2)2 


Si> 


where  gx  €  C\Rd)+,  g2  6  S{Rd),  has  a  unique  strong  solution  on  R+  ,  such  that  U{t),  dU  ( t)/dt ,  jA  U(t) 
are  continuous  functions  from  R+  to  C\Rd),  and  U{t)  €  C\Rd)jr  for  each  t  €  R+  . 
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Recall  that  C{.ftd)+  is  defined  by  (4.41). 

Lemma  A.2  For  each  gi  €  C'(Rd)+  there  exists  a  unique  strong  solution  of  (7 r)  on  [0,  tmai), 
where  tmax  <  00.  Moreover,  iftmax  <  00,  then  lim  ||t/(l)||  =  00. 


Proof  The  proof  was  outlined  by  Perkins  [12]  for  the  case  of  real- valued  functions.  In  our  case 
the  proof  is  completely  the  same  (c.f.  Theorems  (3.13  —  3.15)  from  [12]).  The  more  complicated 
part  of  existence  of  a  solution  without  explosion  will  be  given  here.  That  is,  we  shall  show  that 

tmai  —  OO. 

Before  starting  our  proofs  we  introduce  several  lemmas  from  [12]. 

Lemma  A. 3  Let  Xm  be  the  super- Brownian  motion  given  by  the  unique  solution  of  martingale 
problem  for  (A,£m),  m  €  Mp,  where 

A  =  jexp  {-($,//)},  exp{-($,/i)}(-as+  \{og)2  -  \Ag,p)  :  g  6  T>  (§a)  +  ,  a,o  6  /?}  . 

Then  for  each  rj)  :  [0,oo)  x  Rd  — ►  R,  such  that  ip(s) ,  dip(s)  /ds,  ^Arpis)  are  strongly  continuous 
functions  from  [0;  T]  to  C/(Rd)  we  have 

(A-l)  X?  W))  =  m(m)  +  J* ^  A rP(s)  +  a^(s))  ds  +  Z(  (*) , 

where  Zt  (rj>)  is  the  martingale  with  increasing  process 

(^W),={Am^)2)  ds. 

Lemma  A.4  The  unique  strong  solution  of  the  nonlinear  evolution  equation 

=  -i<T2V(02  +  aV(t)  +  §AV(0 

=  *€P(*A)+ 
satisfies  V(t)  €  Ci{Rd)+  for  each  t  >  0. 

Let 


dV(t) 

dt 

.  V(0) 


(A.2)  U(t)  =  tM0  +  *'tM<) 

be  the  unique  strong  solution  of  ( ir )  on  [0,  tm or),  where  U\  and  U2  are,  respectively,  the  real  and 
imaginary  parts  of  the  solution.  We  can  see  from  (tt)  that  that  U\  and  U2  satisfy 


(A.3) 


=  -ianW  +  ^Om  +  jACf^l) 

+  jont^l)2  -  auUiityi  +  5<*22(ff2)2 
,  W)  =  ffn. 


(A.4) 


Whit) 

dt 


_OllUi(t)l/2(t)  +  t/2(l)qi  +  \AU2{t)  +  Oi2t/l(l)p2 


Ui(0)  =  $12 . 


32 


Let  us  define  for  each  h,  f  €  C(Rd) 

(A.5)  F(h,f)  =  \auh2  -  a\3hf  +  |a22/2. 

It  is  easy  to  check  that 

(A.6)  F(hj){x)>0,  Vxefl*. 

Lemma  A.5  Under  the  definitions  (A. 3),  (A.J)  Ux(t)  €  Ci(Rd)+,  for  all  t  ^  Imax- 


Proof  We  shall  prove  the  result  by  contradiction.  Assume  that  there  exists  a  t  <  tmax,  such  that 
infx£fld  U(t,x)  <  0.  Define 


t'  =  inf  (f  <  tmaz  :  inf  f/(<,x)<ol. 
I  xeRd  ) 


The  solution  U(-)  of  (ir)  is  continuous  on  [0,  /miI),  and  g\  6  C\Rd)+  .  Thus  <*  >  0,  and,  for  each 

s  <  r, 


(A.8)  Ux(s)eC,(Rd)+  . 

Define  xp(s)  =  U\(t‘  -  s).  Obviously  we  have 

(hHs)  =  _dU1(u) 

ds  du  u-t’-s 

Now  let  Xm  be  the  super- Brownian  motion  which  is  the  unique  solution  of  the  martingale  problem 
for  (A*,6m),  m  €  Mf,  where 

A*  =  | exp  {-  (<7in/*)}»  exp  {-  (gn,n)}  (-171511  +  5^11  (Sn)2  -  %&9n,v)  :  9n  G  V  (§a)  +  |  . 

By  Lemma  A.3,  and  equations  (A.3),  (A.5)  we  obtain,  for  all  t  <  t*, 

X?(Ut  (<"-<)) 


=  m(£f,(<*))  +  J Cx?  UxQT  -  s)  +  r hUx(f  -  s)j  ds  +  Z,  (0) 

=  m(Ut(r))  +  £x?  (auUx{t)2)  ds  -  J*X?  ( F(U2(T  -  s),g3))  ds  +  Z,  (tl>) , 
where  Zt(ip)  is  the  martingale  on  [0,t*|  with  increasing  process 

(Zmt  =  j'anXT  (Ui(f  -  s)2)  ds. 

Rearanging,  we  obtain 

-xr  (Ut(f  -  t))  -  fx?  (F(t/a(t*  -  s),S2»  ds  =  -m(U\(tm))  -  Zt  (tf)  -  \  ( Z  (*)), 
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and 


(A-9)  exp  {-*r  (Ui(f  -  0)  -  J‘x?  ( F(U2(t'  -  s),g2))  ds} 

=  exp  {-m(f/,(t*))}  exp  {-Zt  (tf)  -  \  (Z (tf)),}. 

Note  that  exp  {-Zt  (V>)  ~\{Z  W)*}  is>  at  least,  a  local  martingale  on  [0,**],  but  that  the  left  hand 
side  of  expression  (A.9)  is,  by  (A.6),  (A.8),  bounded  and  so  is,  in  fact,  a  martingale.  Setting  t  =  t‘ 
we  obtain 

(A.10)  E  |exp  j-X?  (gu)  -  jf‘  X?  ( F(U2(f  -  s},g2 ))  ds  J  =  exp  {-m(t/,(t*))}. 

On  the  other  hand,  by  (A.6), 

(A.ll)  E^xP{-Xp(9ll)- £  X?(F(U2(f-s),g2))ds}  <  £[exp  {-A?  (gu)  ds}} . 

Define  V(t)  to  be  the  unique  strong  solution  of 

< 

V(0)  =  9n  e£(£A)+. 

By  Lemma  A.4  we  have  that  V(t)  €  Ci(Rd)+  for  each  t  >  0.  By  standard  arguments  from  the 
theory  of  superprocess  we  find 

(A.12)  £J[exp  {-*?(*,,)}]  =  exp{-m(V(f))}. 

Combining  (A. 10),  (A.ll),  (A.12)  we  obtain  that 

exp{-m(i/i(C))}  <  exp{-m(K(t*))}. 

The  above  expression  holds  for  all  m  €  MF.  Therefore,  setting  m  =  6X,  we  obtain 

exp{-f/i(C,*)}  <  exp{-K(t*,x)},  Vx  €  Rd 

and  hence 


infui{t%x)  >  inf  V(f,x). 

ten*  x  eR* 


But  infx€R<  V(t*,x)  >  0,  so  that  infx€/}*  Ui(tm,x )  >  0,  which  is  a  contradiction. 

Lemma  A.6  There  exist  constants  C\ ,  Cj  >  0  such  that 
(A.13)  ll^(0ll<C1exp{C2t},  V t  <  tmax 

and  hence,  if  tmax  <  oo,  then 

(A.14)  limsup  ||(/2(t)||  <  <7,  exp  {Catmar}. 
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Proof  Consider  two  cases 
(1)  j/i  >  0 
Define 


(A. 15) 

By  (A.4)  we  obtain 

dU2(t) 

dt 


u2(t)  =  U2{t)exp{-mt},  t<  tmax  • 

(-anUi(t)U2(t)  +  U2(t)m  +  \AU2(t)  +  a12Ui(t)g2)e\p  {-^t) 
-  rnU2(t)e\p 


=  -aulh(t)U2(J)  +  \AU2(t)  +  al2Ut(t)g2exp  {-ijiO 


=  Ui{t)(-anU2(t)  +  al2g2exp{-r)it}S)  +  \AU2(t),  t  <  tr 

(lM|, 


Define 

C\  =2  max 

Assume  that  there  exists  t0  <  tmax,  such  that 


a22g2 


)■ 


sup  \\u2(t)\\  >  Ci. 


Choose  <*  such  that 

(a. i6)  \\u2(n\\  =  sup  IMI- 

By  choice  of  Ci,  we  have  that  <*  >  0.  For  all  t  <  tmax,  U2(t)  €  S{Rd),  so  that  there  exists  a 
x*  €  Rd,  for  which 

p2(f,**)|  =  ||t/2(C)|. 

Assume,  without  loss  of  generality,  that  x*  is  the  point  of  maximum  of  U2(t‘)  (in  the  case  of  a 
minimum  the  proof  is  analogous).  By  the  previous  lemma,  U\{t)  €  Ci(Rd)+,  for  all  t  <  tm0x-  By 
the  choice  of  C\  we  obtain  that 

l/i(i*,x*)  (-ani?2(t*,x*)  +  ai2$2 exp  <  0. 


Since  x*  is  the  point  of  maximum  of  U2(tm),  the  positive  maximum  principle  ([5],  p.165)  implies 
that 

iA£/2(t*,x’)  <  0. 

Hence 

dr 
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Recall  that  /*  >  0,  so  that  there  exists  a  i  <  tm  for  which  U2(t,x‘)  >  U2(tm,xm).  But  this  contra¬ 
dicts  (A.  16).  Hence  we  obtain  that 


(<sup  |l72(0||  <  Ci, 


and  by  the  definition  of  U2  we  have 

(A.17)  ||t/2(0l|<C1exp{r?1t},  t 

^  fmar  > 

which  was  what  we  wanted  to  prove. 

(2)  T7,  <  0 

The  proof  is  the  same,  the  only  difference  being  that  we  need  not  introduce  f/2,  as  in  the  previous 
case.  B 


Proof  of  Lemma  A.l  Assume  by  contradiction  that  tmax  <  oo.  It  is  well-known  [11]  that  if 
U(t)  is  a  strong  solution  <u  (tt)  on  [0,  tmax)  then  it  is  also  a  mild  solution  of  the  following  integral 
equation  : 

U(t)  —  S(t)g\  +  S(t  —  s)  +  U(s)t]\  +  iaj2[/(s)g2  +  |a22  (<72)2)  ds , 

where  S(t)  is  the  semigroup  generated  by  -A.  By  (A.2),  (A. 3)  we  obtain 
(A.18)  Ui(t)  =  S(t)gn+ £s{t-s)(-\auUi(s)2  +  Ul(s)m  +  F(U2(s),g2j)  ds 
Recall  that,  by  Lemma  A.5,  we  have  U\{i)  €  C\Rd)+  for  all  t  <  tmax.  Thus 
0  <  Ui(t)  =  S(t)gu  +  S(t  -  s)  (-jOnf/iCs)2  -(-  U\(s)t/i^  ds  +  J  S{t  -  s)F  (U2(s),g2)  ds 

<  S(t)gu  +  [*S(t  -  rn)  ds  +  ['s(t  -  s)F(U2(s),g2)  ds 

Jo  Jo 

<  llffnll  +  Irnl  f*  ll^i(s)||  ds  +  t  sup  ||F  (C72(s),  <72)|| 

Jo  $<t 

Finally,  defining  C  =  tmax sup,<<moi  |jF(t/2(s),<72)||,  (C  <  oo  by  Lemma  A.6)  and  K  =  ||<7„||  +  C 
we  obtain 

(A-l»>  lltMOII  <  tf  +  M  fwUiWWds,  vt  <tmax, 

Jo 

and  by  Gronwall’s  inequality  we  have  that 

lltMOII  <  A'exp{|r?,|t},  V«  <  7max  • 

This  implies  that 

(A.20)  lim sup  HtMOU  <  KBXP{\m\tmax}  <  oo. 
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Combining  the  above  result  with  Lemma  A. 6  we  obtain 


(A.21)  lim  ||l/(f)j|  <  limsup||[/i(t)|j  +  lim  sup  ||t/2(t)||  <  °°> 

max  t|fmax 

which  contradicts  Lemma  A.2.  Hence  tmax  =  oo,  which  finishes  the  proof  of  the  lemma. 
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